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Abstract.  An  extensive  airborne  expendable  bathythermograph  survey  of  the 
South  China  Sea  (SCS)  conducted  in  May  1995  and  historical  data  are  used  to 
analyze  and  infer  the  upper  layer  (300  m)  synoptic  structure  and  general  circulation. 
The  primary  thermal  feature  observed  was  a  central  SCS  warm  pool  surrounded 
by  several  cool  pools.  The  size  of  the  warm  pool  decreased  with  depth  from 
approximately  200.000  km2  at  50  m  depth  to  about  70.000  km2  at  300  m  depth. 
The  maximum  temperature  of  the  warm  pool  was  30°C,  appearing  near  the  surface. 
At  the  depth  of  50  m.  the  temperature  of  the  central  SCS  warm  pool  was  29°C, 
and  the  temperature  of  the  five  surrounding  cool  pools  ranged  from  26°C  to  22°C. 
A  three-dimensional  estimate  of  the  absolute  velocity  field  was  obtained  from  the 
observed  temperature  field  and  a  climatological  salinity  field  using  the  3  spiral 
method.  Striking  circulation  features  were  the  existence  of  dual  anticvdonic  eddies 
in  the  centra]  SCS  warm  pool  and  the  existence  of  cyclonic  eddies  associated  with 
the  cool  pools.  In  the  upper  layer  the  tangential  velocity  of  the  dual  central  SCS 
ant  icy  clonic  warm-core  eddies  is  around  30-40  cm/s  and  that  of  the  five  cyclonic 
cool-core  eddies  varies  from  10  cm/s  to  40  cm/s.  The  tangential  velocity  of  ail  the 
eddies  decreased  with  depth.  At  300  m  depth,  it  became  less  than  5  cm/s  for  all 
the  eddies. 


1.  Introduction 

The  South  China  Sea  (SCS)  has  a  bottom  topography 
(Figure  1)  that  makes  it  a  unique  semi-enclosed  ocean 
basin  that  is  temporally  forced  by  a  pronounced  mon¬ 
soon  wind.  Extended  continental  shelves  (less  than  100 
m  deep)  exist  along  the  north  boundary  and  across  the 
southwest  portion  of  the  basin,  while  steep  slopes  with 
almost  no  shelf  are  found  along  the  eastern  boundary. 
The  deepest  water  is  confined  to  an  oblate  bowl  oriented 
SW-NE,  centered  around  1 3° N .  The  maximum  depth  is 
around  4500  m. 

The  SOS's  physical  oceanography  has  been  investi¬ 
gated  for  several  decades.  Studies  show  that  eddies 
in  SCS  are  predominantly  cyclonic  in  winter  and  an- 
ticvclonic  in  summer,  with  sizes  ranging  from  small  to 
mesoscale,  and  both  cold  and  warm  eddies  exist  in  SCS 
[Huang  el  al.,  1994].  A  survey  by  Wyrtki  [1961]  revealed 
complex  temporal  and  spatial  features  of  the  surface 
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currents  in  both  the  SCS  and  the  surrounding  waters. 
By  April,  the  northeast  monsoon  has  decayed,  and  the 
main  SCS  current  consists  of  two  large  cyclonic  eddies  in 
the  western  SCS,  one  in  the  north  and  the  other  in  the 
south  (Figure  2a).  With  the  beginning  of  the  south¬ 
west  monsoon  in  May  and  June,  a  wind  drift  current 
is  formed.  Off  the  coast  of  Vietnam,  westward  inten¬ 
sification  of  this  current  is  clearly  visible.  The  wide, 
uniform  drift  in  the  northern  SCS  shows  a  deflection 
of  the  current  to  the  right  of  the  wind.  A  large  part 
of  the  water  passes  south  of  Taiwan  and  into  the  root 
of  the  Kuroshio.  The  southern  SCS  is  occupied  by  an 
anticyclonic  eddy  (Figure  2b). 

Seasonal  occurrence  of  SCS  eddies  have  been  reported 
by  several  authors.  Dak  [1956]  and  Uda  and  Nakao 
[1972]  reported  a  cold  eddy  off  the  central  Vietnam 
coast  in  summer.  NHani  [1970]  found  a  cold  eddy  lo¬ 
cated  northwest  of  Luzon  in  summer.  Reports  from 
the  Souih  China  Sea  Institute  of  Oceanology  (SCSIO) 
[1985]  indicate  that  a  warm-core  eddy  appears  in  sum¬ 
mer  and  winter  in  the  central  SCS,  but  in  summer,  it 
is  closer  to  Vietnam  at  the  surface.  Recently,  a  cold- 
core  eddy  was  detected  in  the  central  SCS  during  De- 
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Figure  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the 
South  China  Sea. 


cember  29,  1993.  to  JanuaTv  5.  1994,  from  the  analy¬ 
sis  of  TOPEX/  POSEIDON  data  [Soon^  et  al.  1995]. 
Chu  et  at.  [1997a]  and  Chu  and  Chang  [1997]  identi¬ 
fied  the  existence  of  a  central  SCS  surface  warm-core 
eddy  in  mid-Mav  from  a  historical  data  set:  the  U.S. 
Navy's  Master  Observational  Oceanographic  Data  Set 
(MOODS).  From  the  composite  analysis  on  the  U.S.  Na¬ 
tional  Centers  for  Environmental  Prediction  (NCEP) 
monthly  sea  surface  temperature  (SST)  fields  (1982- 
1994).  Chu  et  al.  [1997b]  found  that  during  the  spring- 
to-summer  monsoon  transition  (March-May)  a  warm 
anomaly  (greater  than  1 .8°C)  is  formed  in  the  centra! 
SCS  at  1120-119°30’E,  15®-19#30*N. 

All  of  these  analyses  were  based  on  the  sea  surface 
data.  What  is  the  subsurface  thermal  structure  and  cir¬ 
culation  pattern  of  the  SCS  eddies?  The  Naval  Oceano¬ 
graphic  Office  conducted  an  intensive  airborne  expend¬ 
able  bathythermograph  (AXBT)  survey  between  May 
14-25.  1995.  over  the  majority  of  the  SCS  down  to 
about  300  m  depth.  Figure  3  shows  the  daily  AXBT 
deployment.  This  data  set  provides  something  close  to 
a  "snapshot"  of  the  temperature  in  the  upper  ocean 
in  the  SCS  during  the  transition  time  before  the  on¬ 
set  of  the  monsoon.  Since  only  temperature-measuring 
AXBTs  were  used,  no  salinity  measurements  were  made 


at  the  same  time.  We  will  invert  three-dimensional  eddy 
features  from  analysis  of  this  temperature  data  set  us¬ 
ing  climatological  salinity  data  [Xevtfus,  1984]  in  lieu  of 
in  situ  salinity  data. 

The  outline  of  this  paper  is  as  follows.  A  descrip¬ 
tion  of  the  AXBT  measurement  and  establishment  of 
gridded  data  is  given  in  section  2.  A  depiction  of  the 
3  spiral  method  with  two  necessary  conditions  is  given 
in  section  3.  The  synoptic  three-dimensional  thermal 
structure  and  inverted  velocity  field  is  discussed  in  sec¬ 
tion  4.  In  section  5  we  present  our  conclusions. 

2.  AXBT  Measurements 

Most  of  the  376  AXBTs  were  deployed  at  six  intervals 
over  a  12  day  period  from  May  14-25.  1995  (Figure  3). 
The  majority  of  the  AXBTs  were  nominally  capable  of 
reaching  a  depth  of  360-400  m.  The  ensemble  of  tem¬ 
perature  profiles  (Figure  4a)  and  the  mean  profile  with 
an  envelope  of  a  standard  deviation  (Figure  4b)  show 
the  existence  of  a  mixed  layer  with  depths  ranging  from 
20  to  60  m  and  a  thermocline  with  a  vertical  temper¬ 
ature  gradient  of  6°-7®C  per  100  m  below  the  mixed 
layer. 

The  AXBT  observations  were  mapped  to  a  regular 
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Figure  2.  A  drastic  change  of  the  surface  circulation  from  (a)  April  to  (b) 
June  [after  Wyrtki,  1961]. 


grid  at  12  levels  (Table  1)  using  a  two-scale  optimal  in¬ 
terpolation  (01)  scheme  [Gandtn,  1965;  Lozano  et  al., 
1996].  The  large-scale  01  was  used  to  estimate  the 
background  mean  with  a  decorrelation  scale  of  450  km. 
The  mesoscale  01  was  used  to  map  the  observational 
anomaly  from  the  background  mean  field  into  a  reg¬ 
ular  grid  with  a  spatial  decorrelation  scale  of  75  km 
and  a  temporal  decorrelation  scale  of  10  days.  The 
spatial  and  temporal  decorrelation  scales  associated  to 
the  mean  field  was  estimated  from  the  covariance  ma¬ 
trix  obtained  from  MOODS  for  the  SCS  [see  Chv  et 


al..  1997c].  The  large-scale  mean  was  computed  using 
the  temperature  and  salinity  profiles.  The  large-scale 
pattern  for  the  surface  (not  shown)  resembles  the  pat¬ 
tern  previously  estimated  by  Chu  el  at.  [1997c],  The 
decorrelation  length  scale  for  the  fluctuations  from  the 
mean  was  estimated  from  the  AXBT  data.  The  results 
shown  below  are  relatively  insensitive  to  the  length  scale 
in  the  range  of  60-90  km.  The  decorrelation  timescale 
was  selected  to  ensure  synopticity.  The  SST  horizontal 
scales  seem  somewhat  larger  than  temperatures  in  the 
seasonal  thermocline.  In  this  study,  we  choose  horizon- 
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Figure  3  Deployment  pattern  of  AXBT  survey  during  May  14-25.  1995. 


tal  scales  to  be  uniform  in  the  vertical  for  convenience 
and  for  lack  of  sufficient  data  to  discriminate  these  dif¬ 
ferences  properly. 

3.  3  Spiral  Method 

3.1.  Geostrophic  Velocity  Inverted  from 
Density  Field 

Often  the  thermal  wind  relation  is  employed  to  cal¬ 
culate  geostrophic  velocity  from  hydrographic  data 


9 

U  =  Uo+  — 
fpo 

/> 

(1) 

9 

V  =  Vo  -  7 — 
JPo  . 

jO 

(2) 

where  (u,  v)  and  (uq,  uq)  are  the  geostrophic  velocity  at 


any  depth  z  and  at  a  reference  depth  zo,  p  is  the  in 
situ  water  density,  po  is  the  characteristic  value  of  the 
density.  /  =2QsinA,  is  the  Coriolis  parameter,  where 
fi  is  the  Earth’s  rotation  rate  and  A  is  the  latitude. 
The  hydrographic  data  determine  the  geostrophic  ver¬ 
tical  shear.  The  reference  velocity  (u0,t'o)  remains  un¬ 
determined.  On  the  basis  of  the  geostrophic  and  hydro¬ 
static  balances  and  mass  conservation  equations,  sev¬ 
eral  inverse  techniques  have  been  developed  to  compute 
the  geostrophic  velocity  at  the  reference  level  (u0,t>o)- 
These  include  the  3  spiral  method  [5/omme/  and  Schott, 
1977;  Schott  and  Siommel,  1978],  the  box  inverse  method 
[Wunsch,  1978],  and  the  Bernoulli  method  [Killworth, 
1986].  Wunsch  [1994]  proposed  a  generalized  3  spiral 
and  geostrophic  balance  inverse  method  based  on  a  non¬ 
linear  optimization  on  Needler’s  P  equation  [Needier, 
1967].  (See  the  review  paper  on  the  3  spiral  method  by 
Others  et  at.  [1985].)  Davis  [1978]  pointed  out  that  the 
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(a) 


(b) 


Figure  4.  Temperature  profiles  of  AXBT  survey  during  May  14-25,  1995. 
(a)  Ensemble  of  profiles  and  (b)  ihe  mean  profile  with  an  envelope  of  a  st  andard 
deviation.  Both  show  the  existence  of  a  mixed  layer  with  depths  ranging  from  20 
to  60  m  and  a  thermocline  with  a  vertical  temperature  gradient  of  6°-7°C/100 
m  below  the  mixed  layer. 


3  spiral  method  and  the  box  inverse  method,  no  mat¬ 
ter  how  different  in  appearance,  are  based  on  the  same 
dynamical  assumptions,  and  the  existence  of  the  3  spi¬ 
ral  thus  becomes  a  required  condition  for  these  inverse 
methods. 

Recently,  nongeostrop  hie  inverse  methods  were  also 
developed.  Bogden  et  al.  [1993]  estimated  the  time- 
averaged  velocity  field  in  the  North  Atlantic  from  ob¬ 
servations  of  density,  wind  stress,  and  bottom  topog¬ 
raphy  based  on  the  assumption  of  minimum  mixing  of 
potential  density  at  middepth.  Tzipcrman  et  al.  [1992] 
applied  the  adjoint  method  to  investigate  the  steady 
state  oceanic  general  circulation  using  a  realistic  numer¬ 
ical  oceanic  general  circulation  model  and  hydrographic 
data. 

As  pointed  out  by  Wvnsch  and  Grant  [1982],  in  deter¬ 
mining  large-scale  circulation  from  hydrographic  data, 
we  can  be  reasonably  confident  of  the  assumptions  of 
geostrophic  balance,  mass  conservation,  and  no  major 
cross-isopycnal  mixing  (except  for  water  masses  in  con¬ 
tact  with  the  atmosphere).  Under  these  conditions,  the 
density  of  each  fluid  element  would  be  conserved,  which 
mathematically  is  given  by 

V  •  Vp  =  0  (3) 

where  p  is  the  potential  density  and  V  =  (u.v.  w)  is 
the  geostrophic  velocity.  The  conservation  of  potential 
vorticity  equation  [Pedlosky.  1986]  can  be  obtained  by 
differentiating  (3)  with  respect  to  z,  using  geostrophic 
and  hydrostatic  balances,  and  including  the  latitudinal 
variation  of  the  Coriolis  parameter  to  give 

V  •  V?  =  0  (4) 


where  q  =  fdp/dz.  Use  of  fdp/dz  may  induce  a  small 
but  systematic  error  into  estimation  of  potential  vor¬ 
ticity  [G.T.  Needier,  unpublished  manuscript.  1986]. 
Equations  (3)  and  (4)  imply  that  the  velocity  V  is  par¬ 
allel  to  Vq  x  Vp. 

3.2.  Necessary  Conditions  for  Validity  of  the  .3 
Spiral  Method 

Stommel  and  Schott  [1977]  pointed  out  that  the  three- 
dimensional  velocity  field  cannot  be  determined  from 
the  density  field  alone  when  the  q  and  p  surfaces  coin¬ 
cide.  The  first  necessary  condition  for  the  validity  of 
this  inverse  method  thus  is 

1.  Necessary  Condition  1:  The  p  surface  is  not  par¬ 
allel  to  the  q  surface,  which  mathematically  requires 

Vp  x  Vq  £  0  (5) 

Stommel  and  Schott  [1977]  further  pointed  out  that 
the  three-dimensional  velocity  field  cannot  be  deter¬ 
mined  from  the  p  field  alone  if  the  horizontal  velocity 
does  not  rotate  with  depth  {3  spiral).  The  existence  of 
the  3  spiral  is  the  second  necessary  condition. 

2.  Necessary  Condition  2:  The  velocity  («,  v)  should 

execute  a  3  spiral,  which  mathematically  requires  that 
for  at  least  two  depths,  z  =  s*.  and  z  =  zm,  with  hori¬ 
zontal  velocities  «<*>]  and  [u(m).t;(m)]. 


If  we  cannot  find  levels  z*  and  zm  such  that  (6)  is  satis¬ 
fied.  the  inverse  method  will  fail  to  get  velocity  in  that 
water  column. 
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Table  1.  Vertical  Levels  of  the  Interpolated  Tempera¬ 
ture  Field 


Level 

Depth  (m) 

1 

0 

2 

10 

3 

20 

4 

30 

5 

50 

6 

75 

7 

100 

s 

125 

9 

150 

10 

200 

11 

250 

12 

300 

3.3.  Implementation  of  the  Necessary 
Conditions 

Consider  the  unit  vector  P  [Chu.  1995]  defined  by 


Vp  x  Vq 
|Vp  x  Vq\ 


(7) 


The  existence  of  this  unit  vector  implies  the  satisfaction 
of  necessary  condition  1. 

The  velocity,  V  =  (u,v,w).  parallels  the  unit  vector 

P. 


V=r(*,y,*)P  (8) 

where  r  is  the  proportionality.  Applying  the  thermal 
wind  relation  at  two  different  depths  zt  and  zm,  as 
shown  in  Figure  5a.  a  set  of  algebraic  equations  for  de¬ 
termining  the  parameter  r  are  obtained 


r(*>pi:*>  _  r(^)p(m)  =  Avkm  (9) 

which  are  two  linear  algebraic  equations  for  and 
r<m:  [rW  =  r(ar,  y,  z,-)] .  Here 


(AUkm'AVkm)~  fpa  Lm(di  dx)d~  (10) 
The  determinant  of  (9) 


pjT'  Px”^ 

p(*)  p(m) 

rv  rv 


=  sin(orjtm) 


(11) 


(k) 

is  the  sine  of  the  vertical  turning  angle  between  Ph  and 

P'hm>  (Figure  5b);  that  is.  o*m  is  the  3  spiral  turning 
angle  between  the  two  levels  r*  and  zm.  The  algebraic 
equations  in  (9)  have  definite  solutions  for  r(  *1 


if 


= 


P‘m) 
Artm  Pjm} 

sin(a*m) 


(12) 


a*m^0  (13) 

As  soon  as  r<*>(*=1.2,  ...)  is  obtained,  the  absolute 
velocity  field  can  be  calculated  by  (8).  Before  using  the 
inverse  method,  we  should  check  for  the  two  necessary 
conditions  (5)  arid  (6)  or.  equivalently,  (7)  and  (13). 

For  water  columns  satisfying  the  two  necessary  con¬ 
ditions.  we  may  use  (13)  to  compute  r(i'  for  the  level 
There  are  A‘-l  sets  (m  =  1,  2,  Jt-1,  Jt+.l,  ...  .V)  of 
equations  (9)  for  calculating  r[t\  Here  N  is  the  total 


Z 


Figure  5.  Vertical  turning  of  the  P  vector,  (a)  P  vector  at  two  different 
levels  and  (b)  turning  angle  between  two  levels. 
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Figure  6.  Horizontal  temperature  fields  at  different  depths:  (a)  0  m,  (b) 
50  tn,  (c)  100  m,  (d)  200  m,  and  (e)  300  m.  The  central  SCS  warm  pool  was 
surrounded  by  Dongsha  (DS),  Hainan  (HN),  South- Vietnam  (SV),  and  Liyue 
(LY)  cool  pools  and  Xisha  (XS)  warm  pool. 


number  of  vertical  levels  of  the  water  column.  All  the 
.V-l  sets  of  equations  are  compatible  under  the  ther¬ 
mal  wind  constraint  and  should  provide  the  same  so¬ 
lution.  However,  because  of  errors  in  measurements 
(instrumentation  errors)  and  computations  (truncation 
errors),  the  parameters  r(k)  may  vary  with  m.  A  least 
squares  error  algorithm  (appendix  A)  is  used  to  mini¬ 
mize  the  error.  For  further  details  and  validation  of  the 
algorithm  presented  here  sec  Chu  ct  al.  [1998].  Inter¬ 
ested  readers  can  obtain  the  software  by  contacting  Mr. 
Chenwu  Fan  at  fanSnps.navy.mil  or  visiting  our  Web 
site  at  http://web.nps.navy.mil/'chu. 


4.  Results 

4.1.  Thermal  Field 

4.1.1.  Horizontal  structures.  Figures  6a-6e  are 
horizontal  depictions  of  temperature  at  0.  50.  100,  200. 
and  300  rn  depths,  respectively.  The  contour  interval 
is  0.5°C.  In  them,  we  see  that  warmer  water  is  situ¬ 
ated  at  the  centra)  SCS  with  surrounding  cooler  water. 
At.  the  50  m  depth,  the  28°C  isotherm  separates  the 
warm  central  SCS  water  from  the  surrounding  cooler 
water.  The  maximum  temperature  reaches  29° C.  This 
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Table  2.  Locations,  Typical  Temperatures,  and  Tangential  Velocity  of  the  SCS  Warm-Core  and  Cool-Core 
Eddies  Appearing  During  May  14-25.  1995 


Depth  (m) 

Central  SCS  (W) 

Xisha  (W) 

Dongsha  (C) 

Hainan  (C) 

South  Vietnam  (C)  Liyue  (C) 

0 

1120-117®E 

110®-1 12®E 

1 16°-118°20'E 

110®-112®E 

108°-112°£ 

1 15°- 1 18°  F. 

I2°-1S°.N 

14®-16®N 

18®-20°20'N 

18°-20®N 

10°-14°N 

12°-!4°N 

29.5°-30°C 

29.5°C 

27®C 

28®C 

28.5®C 

29.5®-30®C 

40  cm/s 

20  cm/s 

40  cm/s 

10  cm/s 

20  cm/s 

30  cm/s 

50 

112®-117°E 

110®-112®£ 

116°-118°20'E 

I10®-112®E 

108®-112®E 

115°-118°E 

12°-18°N 

14®-16°N 

18®-20°20'N 

16®-18®N 

10“-14°N 

12°-14°N 

29°C 

24®  C 

22°C 

23®  C 

23°  C 

26°C 

30  cm/s 

10  cm/s 

25  cm/s 

8  cm/s 

10  cm/s 

30  cm/s 

100 

112°-1 17®E 

110®-112®E 

1 16°-1 18®20'E. 

109®-111°30'E 

110®-112®E 

1 15®-118®E 

12®-18®N 

14°-13®N 

1S®-20®20'N 

16®40'-19°N 

12®-14®N 

12®-14°N 

23®C 

21°C 

18®C 

18®C 

20®  C 

19®C 

25  cm/s 

8  cm/s 

20  cm/s 

5  cm/s 

10  cm/s 

20  cm/s 

200 

112®-117®E 

110®-112®E 

116®-118®20'E 

109®-113°E 

1 14°-117®E 

12®-18®N 

14®-15®N 

]8®-‘20®20'N 

16®-18°N 

11®-14°N 

16®C 

15“C 

13®C 

13.5®C 

14°C 

10  cm/s 

4  cm/s 

10  cm/s 

4  cm/s 

9  cm/s 

300 

115®-117®E 

... 

1 16®-118°20'E 

1 10®-113®E 

114°-116®E 

16°-18®N 

.  .  .. 

18°-20®20'N 

16°-18°N 

t.r 

I0#-14*N 

12®C 

11®C 

11®C 

1 1.5*C 

2  cm/s 

.  .  . 

3  cm/s 

1  cm/s 

3  cm/s 

central  SCS  warm  pool  extended  to  the  300  m  depth 
and  was  surrounded  by  Dongsha  (DS),  Hainan  (HN), 
South- Vietnam  (SV),  and  Liyue  (LY)  cool  pools  and 
Xisha  (XS)  warm  pool.  The  location  and  thermal  fea¬ 
tures  of  these  warm  and  cool  pools  are  listed  In  Table 
2. 

4.1.2.  Vertical  structures.  Four  zonal  cross  sec¬ 
tions  (19®  N .  1 7®  N .  13®  N .  and  7®  N)  of  temperature  show 
the  vertical  structure  of  the  mixed  layer  and  the  thermo- 
cline  as  well  as  the  warm  and  cool  pools  (Figure  7).  The 
mixed  layer  was  not  evident  at  the  I9°N  cross  section 
but  develops  as  latitude  decreases.  For  the  17®  N  and 
7®N  cross  sections,  the  mixed  layer  shallows  toward  the 
east.  In  the  13° N  cross  section,  the  mixed  layer  deepens 
at  114®E.  The  trough  and  ridge  of  the  isotherms  may 
be  considered  as  centers  of  warm  and  cool  pools. 

The  19®N  cross  section  clearly  shows  the  existence  of 
DS  cool  pool  near  117°E  from  the  uplifting  (ridge)  of 
isotherms.  The  ridge  is  located  at  117®E  at  the  surface, 
eastward  shifted  with  depth,  and  kept  at  117°E  below 
the  100  m  depth.  The  longitudinal  span  of  DS  cool  pool 
is  around  2®. 

The  two  cross  sections  at  17°N  and  13® N  show  the 
occurrence  of  a  central  SCS  warm  pool  from  the  down¬ 
ward  bending  (trough)  of  isotherms:  1I6®E  at  17®N  and 
114°E  at  13® N.  The  longitudinal  span  of  central  SCS 
warm  pool  is  around  4®  for  both  latitudes.  The  warm 
water  is  nearly  isolated  from  the  surrounding  cooler  wa¬ 
ter.  This  may  suggest  that  local  forcing  might  be  im¬ 
portant  for  the  central  SCS  warm  pool  formation. 

The  7®N  cross  section  shows  the  eastward  uplift  of 
the  thermocline  depth  from  50  m  at  110®E  to  25  m 
at  115®E.  A  surface  warm  pool  with  SST  higher  than 


30® C'  was  found  between  112®-115®E,  and  a  midlevel 
(100  rn)  cool  pool  was  detected  between  110°-112®E.  In 
the  eastern  part  (1 13®-1 15®E),  an  isolated  warm  pool 
appeared  below  the  175  m  depth. 

Three  latitudinal  cross  sections  (111°E,  114®E,  and 
117®E)  of  temperature  also  show  the  vertical  structure 
of  the  mixed  layer  and  the  thermocline  as  well  as  the 
warm  and  cool  pools  (Figure  8).  The  mixed  layer  had 
a  shallow-deep-shallow  pattern  in  the  latitudinal  direc¬ 
tion.  Taking  the  114®E  cross  section  as  an  example, 
the  mixed  layer  had  a  shallow  depth  (about  10  m)  near 
5®N.  Its  depth  increased  with  latitude  to  about  60  m  at 
13°30’N  and  decreased  with  latitude  to  almost  zero  (no 
mixed  layer)  near  20® N. 

The  ‘‘trough”  and  “ridge”  of  the  isotherms  may  be 
considered  as  centers  of  warm  and  cool  pools.  The 
three  cross  sections  indicated  the  occurrence  of  cen¬ 
tral  SCS  warm  eddy  evident,  from  the  downward  bend¬ 
ing  (trough)  of  isotherms:  9®N  at  the  111®E,  13®N  at 
1 14®E.  and  15®N  and  17®N  at  117®E. 

4.2.  The  SCS  Circulation 

4.2.1.  T-S  and  T-S  relations.  The  3  spiral 
method  was  used  to  determine  the  velocity  field  from 
hydrographic  data.  Since  there  were  no  salinity  (5)  ob¬ 
servations  in  SCS  during  May  14-25,  1995.  May  clima¬ 
tological  salinity  [Levitus.  1984]  was  interpolated  to  the 
AXBT  stations  and  those  salinity  values  C5)  were  used 
for  the  P  vector  computation.  To  validate  this  tech¬ 
nique.  we  plotted  two  T-S  diagrams:  One  diagram  (Fig¬ 
ure  9a)  is  the  plot  of  3478  T.  S  profiles  from  MOODS 
in  May  between  1933-1986  and  the  other  diagram  (Fig- 
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Figure  7,  Temperature  distribution  at  several  zonal 
cross  sections:  (ai  19SN,  (b)  17°N,  (c)  13°N,  and  (d)  7®N. 


ure  9b)  is  the  plot  of  376  T  (AXBT)  and  S  profiles. 
The  spatial  and  temporal  distributions  of  the  MOODS 
data  are  illustrated  in  Figures  10a  and  10b.  Comparison 
between  the  two  diagrams  shows  that  the  T-S  charac¬ 
teristics  remain  stable.  However,  we  should  be  aware 
of  the  shortcomings  in  using  the  climatological  salinity 
values  (5).  If  the  temperature  field  is  relatively  homo¬ 
geneous  in  space,  the  density  gradient  depends  mostly 
on  the  salinity  gradient.  Under  that  circumstance,  the 
use  of  the  climatological  salinity  values  may  bring  large 
errors  in  inverting  the  velocity  field. 

4.2.2.  Horizontal  velocity  fields.  The  available 
historical  data  [Wyrtkt.  1961]  show'  two  totally  different 
circulation  patterns  for  the  months  preceding  and  fol¬ 
lowing  the  AXBT  survey  in  the  central  SCS:  a  cyclonic 
gyre  in  April  (Figure  2a)  and  an  anticyclonic  gyre  in 
June  (Figure  2b).  We  may  treat  the  May  circulation  as 
the  transition  between  the  two  patterns. 

Our  inverted  velocity  data  set  (Figure  11)  indicates 
the  existence  of  a  central  SCS  anticyclonic  gyre  sur¬ 


rounded  by  several  cyclonic  and  anticyclonic  eddies. 
Generally,  the  anticyclonic  (cyclonic)  eddies  are  asso¬ 
ciated  with  the  warm  (cool)  pools.  The  kinetic  features 
of  these  eddies  are  listed  in  Table  2.  Furthermore,  the 
magnitude  of  inverted  velocities  decreases  with  depth 
and  is  quite  small  at  300  m.  This  implies  that  a  di¬ 
agnostic  met  hod  with  a  level  of  no  motion  at  or  below 
about  a  300  m  depth  could  have  been  used  to  give  re¬ 
sults  close  to  the  3  spiral  method. 

The  observed  three-dimensional  temperature  and  in¬ 
verted  velocity  fields  indicate  the  existence  of  a  central 
SCS  warm-core  anticyclonic  eddy  with  surrounding  DS, 
HN.  SV.  and  LY  cool-core  cyclonic  eddies  and  XS  warm- 
core  anticyclonic  eddy  (Figure  12). 

4.2.3.  Zonal  cross  sections  of  v  component. 
Four  zonal  cross  sections  (19°N,  17°N.  13°N.  and  7°N) 
of  v  velocity  show  the  vertical  eddy  structure  (Fig¬ 
ure  13).  The  positive  values  indicate  northward  ve¬ 
locity,  and  the  negative  values  refer  to  the  southward 
velocity.  Alternate  positive  and  negative  areas  indicate 
the  occurrence  of  cyclonic  and  anticyclonic  eddies.  At 
each  zonal  cross  section,  a  neighboring  eastern  nega¬ 
tive/western  positive  pattern  refers  to  an  anticyclonic 
eddy,  and,  a  neighboring  western  negative/eastern  pos¬ 
itive  pattern  refers  to  a  cyclonic  eddy. 

The  19°N  cross  section  clearly  shows  the  existence  of 
the  DS  cyclonic  eddy  near  1 1 7° E.  The  diameter  of  the 
DS  eddy  is  around  300  km.  The  depth  of  the  eddy  is 
around  200  m.  This  eddy  reveals  an  asymmetric  fea- 
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Figure  8.  Temperat  ure  distribution  at  several  latitudinal 
cross  sections:  (a)  111°E,  (b)  134°E,  and  (c)  1176E. 
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Figure  9.  T-S  diagrams  for  (a)  the  Navy’s  MOODS  dataset  in  May  between 


1933-1986  and  (b)  the  AXBT  temperature  and  May  climatological  salinity. 


ture.  The  northward  velocity  with  a  maximum  speed 
of  40  cm/s  in  the  eastern  part  is  much  stronger  than  the 
southward  velocity  with  a  maximum  speed  of  10  cm/s 
in  the  western  part. 

The  17eN  cross  section  clearly  shows  the  existence  of 
the  centra]  SCS  anticyclonic  eddy  centered  at  116°E. 
The  zonal  span  of  the  central  SCS  eddy  is  around  500 
km.  The  depth  of  the  eddy  is  around  300  m.  The 
maximum  tangential  velocity  is  around  20  cm/s.  This 


eddy  is  quite  symmetric  with  a  maximum  southward 
velocity  in  the  eastern  part  similar  to  the  northward 
velocity  in  the  western  part. 

The  13°.N  cross  section  clearly  shows  the  existence 
of  a  double-eddy  structure.  Between  112°E  and  1 1 6° E 
there  was  a  strong  anticyclonic  eddy  (the  central  SCS 
eddy)  with  a  maximum  tangential  speed  of  30  cm/s  (ap¬ 
pearing  in  the  upper  50  m)  and  centered  at  114°E.  The 
zonal  span  of  this  eddy  is  around  500  km.  This  anti- 
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Figure  10.  MOODS  data  distribution:  (a)  spatial  and  (b)  temporal. 
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Figure  11.  Absolute  velocity  at  different  depths:  (a)  0  m.  (b)  50  m.  (c)  100 
m.  (d)  200  m.  and  (e)  300  m.  The  central  SCS  anticyclonic  eddy  was  surrounded 
by  DS.  HN.  SV,  and  LY  cyclonic  eddies  and  XS  anticyclonic  eddy. 
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cyclonic  eddy  is  associated  with  the  central  SCS  warm 
pool.  The  depth  of  the  eddy  is  around  250  m.  To 
the  east  of  the  anticyclonic  eddy,  a  cyclonic  eddy  ap¬ 
peared  between  114°E  to  119°E  centered  at  1I6°20’E. 
This  cyclonic  eddy  is  associated  with  the  LY  cool  pool, 
and  therefore  called  the  LY  cool-core  (cyclonic)  eddy. 
The  central  SCS  warm-core  eddy  and  LY  cool-core  eddy 
shared  the  same  southward  branch.  Furthermore,  the 
northward  flow  at  the  eastern  part  of  the  IY  cool-core 
eddy  (maximum  tangential  speed  around  10  cm/s)  is 
much  weaker  than  the  southward  flow  in  the  western 
part. 

The  7®N  cross  section  shows  the  occurrence  of  a  cy¬ 
clonic  eddy  (110M14°E).  This  cyclonic  eddy  is  quite 
asymmetric.  The  southward  flow  (maximum  tangential 
speed  around  20  ern/s)  in  the  western  part  is  stronger 
than  the  northward  flow  (maximum  tangential  speed 
around  10  cm/s)  in  the  eastern  part.  The  depth  of  the 
eddy  is  around  150  m.  To  the  east  of  the  anticyclonic 
eddy,  a  cyclonic  eddy  appeared  between  114°E  to  119°E 
centered  at  116°20'E. 

4.2.4.  Latitudinal  cross  sections  of  u  compo¬ 
nent.  Three  latitudinal  cross  sections  (11PE.  114°E, 
and  117°E)  of  a  velocity  also  show  the  vertical  eddy 
structure  (Figure  14).  The  positive  values  indicate  east¬ 
ward  velocity,  and  the  negative  values  refer  to  west¬ 
ward  velocity.  Alternate  positive  and  negative  areas 
indicate  the  occurrence  of  cyclonic  and  anticyclonic  ed¬ 
dies  At  each  latitudinal  cross  section,  a  neighboring 
southern  negative/northern  positive  pattern  refers  to 
an  anticyclonic  eddy.  However,  a  neighboring  northern 


negative/southern  positive  pattern  refers  to  a  cyclonic 
eddy. 

The  area  with  a  southern  negative/northern  positive 
m  pattern  (anticyclonic  eddy)  coincides  with  the  area  of 
the  downward  bending  isotherms  (warm  pool).  At  the 
1 1 1° E  cross  section,  the  central  SCS  warm-core  eddy  is 
recognized  by  an  isotherm  trough  located  at  9°N  (Fig¬ 
ure  12),  and  an  anticyclonic  eddy  is  identified  by  the 
southern  negative/northern  positive  u  pattern  centered 
at  9®N  (Figure  13).  At  the  114°E  cross  section,  the  cen¬ 
tral  SCS  warm-core  eddy  is  recognized  by  an  isotherm 
trough  located  at  13°N  (Figure  12),  and  an  anticyclonic 
eddy  is  identified  by  the  southern  negative/northern 
positive  «  pattern  centered  at  1 3° N  (Figure  13). 

4.3.  Comparison  With  Previous  Work 

This  snapshot  of  multi-eddv  structure  agrees  quite 
well  with  previous  studies  of  individual  SCS  eddies  ap¬ 
pearing  in  mid-May.  For  example,  the  central  SCS 
warm-core  eddy  was  previously  reported  by  the  South 
China  Sea  Institute  of  Oceanology  [SCSfO.  1985]  and 
identified  as  a  recurring  phenomenon  by  analyzing  his¬ 
torical  data  such  as  MOODS  [Chu  et  al..  1997a]  and 
N'CEP  monthly  SST  [Chu  et  al..  1997b].  The  surround¬ 
ing  smaller  cool-core  and  warm-core  eddies  were  identi¬ 
fied  separately  by  different  authors.  The  DS  cool-core 
eddy  was  identified  by  He  et  al.  [1996]  and  the  S\ 
cool-core  eddy  was  reported  by  Dale  [1956]  and  l  da 
and  Nakao  [1972].  Our  contribution  here  is  to  present  a 
synoptic  multi-eddv  structure  in  mid-May,  the  monsoon 
transition  period. 

4.4.  Mechanisms 

In  searching  for  mechanisms  for  the  formation  of  this 
recurring  seasonal  signal,  a  central  SCS  warm  anomaly, 
Chu  et  al.  [1997b]  found  a  high  correlation  between 
surface  wind  stress  curl  and  SST  anomaly  from  correla¬ 
tion  analysis  of  the  European  Center  for  Medium-range 
Weather  Forecast  (ECMVVF)  analyzed  wind  stress  data 
and  the  N'CEP  SST  data.  The  high  correlation  leads  to 
a  hypothesis  for  the  wind-driven  central  SCS  warm  pool 
during  the  spring-to-summer  monsoon  transition.  The 
bowl-type  bottom  topography  shown  in  Figure  1  pro¬ 
vides  a  favorable  condition  for  the  central  SCS  warm 
pool  formation  in  the  spring-to-summer  monsoon  tran¬ 
sition.  In  spring,  a  surface  anticyclone  usually  appears 
over  the  central  SCS  [Cheang,  1987].  Ekman  down- 
welling  will  occur  in  the  cent  ral  part  of  the  "bowl,’’  and 
mass  balance  will  cause  upwelling  near  the  boundary 
of  the  "bowl."  The  down  welling  prevents  the  deep  cold 
water  from  advecting  upward,  while  the  upwelling  along 
the  boundary  helps  the  deep  cold  water  advect  upward. 
This  makes  the  central  part  of  the  "bowl”  warm  and  the 
side  part  of  the  “bowl”  cool  and  leads  to  the  generation 
of  the  central  SCS  warm  pool  surrounded  by  several  cool 
pools.  Furthermore,  other  mechanisms  such  as  remote 
forcing  and  topographic  influence  could  be  significant. 

5. Conclusions 

1.  Objective  analysis  of  the  AXBT  data  allowed  us 
to  identify  the  occurrence  of  a  warm  pool  in  the  central 
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Figure  12  Illustration  of  South  China  Sea  warm-core  and  cool-core  eddies 
in  May  1995. 


SCS  during  May  14-25,  1995.  The  maximum  SST  of 
this  warm  pool  reached  29.5°C  and  higher.  This  warm 
pool  extended  to  the  SCS  southern  boundary  at  the 
surface  and  became  more  isolated  in  the  central  SCS 
below  the  surface.  Its  size  was  near  200,000  km2  at  50 
m  depth  and  decreased  with  depth.  The  double  anti- 
cyclonic  eddies  derived  from  the  inverse  method  were 
colocated  with  the  central  SCS  warm  pool  and  had  a 
maximum  tangential  velocity  around  40  cm/s. 

2.  The  central  SCS  warm  eddy  was  surrounded  by 
several  cool-core  and  warm-core  eddies  with  a  radius 
of  around  100  km  and  maximum  tangential  velocities 
ranging  from  10  cm/s  to  40  cm/s.  The  cool-core  eddies 
were  cyclonic  and  warm-core  eddies  were  anticyclonic. 
The  inverted  velocity  field  shows  the  existence  of  XS 
warm-core  anticyclonic  eddy  with  a  maximum  tangen¬ 
tial  velocity  of  20  cm/s,  strongest  between  the  surface 


and  100  m  depth.  DS  cool-core  cyclonic  eddy  with  a 
maximum  tangential  velocity  of  40  cm/s,  strongest  be- 
tween  the  surface  and  200  ni  depth,  SV  cool-core  cy¬ 
clonic  eddy  between  surface  and  100  m  depth  with  a 
maximum  tangential  velocity  of  20  cm/s  at  surface  and 
decreasing  with  depth,  LY  cool-core  cyclonic  eddy  with 
a  maximum  tangential  velocity  of  30  crn/s.  strongest 
between  the  surface  and  50  m  depth,  and  HN  cool-core 
cyclonic  eddy  with  a  maximum  tangential  velocity  of  10 
cm/s.  the  weakest  of  the  cool-core  eddies. 

3.  The  central  SCS  warm-core  eddy  identified  in 
the  May  1995  AXBT  measurements  is  a  recurring  phe¬ 
nomenon.  The  results  obtained  in  this  study  agree  well 
with  the  historical  data  analyses  such  as  the  MOODS 
and  NCEP  monthly  sea  surface  data  set. 

4.  The  wind-driven  mechanism  for  the  formation  of  a 
central  SCS  warm-core  eddy  with  surrounding  cool-core 


Depth  Depih  Depth  Cteplh 
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level  m.  If  we  assume  that  the  P  vector  exists  at  the 
level  m  and  that  the  velocity  vector  parallels  the  vector 
P',n\  an  error  can  be  easily  defined  by  (Figure  15) 


where  V<m)  =  |V(m>|  and  (P(m>.  V("»>)  is  the  angle 
between  Pfm)  and  V(rn>.  The  total  error  of  the  water 
column  velocity  caused  by  the  uncertainty  of  (u(i\  t4iy) 
is 


where  hm is  the  thickness  of  the  mth  layer  (see  Table  1). 

We  determine  the  velocity  such  that  the 

total  error  E  becomes  minimum;  that  is. 


which  leads  to  a  2x2  algebraic  equations  for  determiu 
ing 

4-  =  F\ 

,421  4  ,422^^  —  E ? 


Loogtituao  (H) 

Figure  13  Distribution  of  velocity  e  component  at  sev 
oral  zonal  cross  sections:  (a)  19°N.  (b)  17®N,  (c)  13°N,  ant 


eddies  during  the  spring-to-summer  transition  was  con¬ 
firmed  by  a  numerical  study.  A  manuscript  on  this  mod¬ 
eling  study  was  submitted  to  the  Journal  of  Oceanog¬ 
raphy  for  publication. 


Appendix:  Optimization  Scheme 

If  the  absolute  velocity  (u(fc),  is  known,  we  may 
use  the  thermal  wind  relation  (10)  to  obtain  the  abso¬ 
lute  velocity  at  any  level  rn. 


dpL^/dz 

The  computed  [uim).  t/m).  u-(m^}  may  not  be  in  the 
same  direction  as  the  P  vector  f Pi-"’ '•  Pirr'  .  P-m/|at  the 


Figure  14.  Distribution  of  velocity  u  component  at  sev¬ 
eral  latitudinal  cross  sections:  (a)  111®E,  (b)  114°E.  and  (c) 
117°£. 
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Figure  15.  Error  caused  by  the  uncertainty  of  the  veloc¬ 
ity  at  the  level  k. 


where 

-4n  =  £m  a(£\  A12  =  Zrn*[? 

••I12  =  A21,  -422  =  Em«22) 

=  ~  Em(all 

and 


a^  =  (P^)^, 


a[™}  = 

4T)  =  (PiTn))2/<m: 

dp'"" /dr  \2  ,  dp^’/Oy  v2j  ,  n 

dpi"'/di  >  J  nm 


Jim  I  ^ 


1  4-  I^LU^L) 2  j_  ( 
1  ’  ''dp^i/dz  >  '  V 


The  absolute  velocity  at  the  level  z  =  z*  can  be  com¬ 
puted  hy 


uft)_  ^1-^22  ~  ^2-412  fyAn  -  FiAn 

-^n-^22  —  -4l2^21  A\\A^Z  —  -•ll2-^2t 

Thus  the  second  necessary  condition  becomes  that  if 

-d  15 -422  —  -4|2-42l  =  0  (14) 

the  inverse  method  fails. 
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